ZU-TH 3/95 
February 1995 
|hep-th/ 9502015 



The number of sphaleron instabilities of the 
Bartnik-McKinnon solitons and non-Abelian 

black holes 

M.S. Volkov, [] O. Brodbeck, G. Lavrelashvili, f] and N. Straumann 

Institut fiir Theoretische Physik der Universitat Ziirich-Irchel, Winterthurerstrasse 190, 

CH-8057 Zurich, 
Switzerland 



Abstract 



It is proven that there are precisely n odd-parity sphaleron-like unstable modes of 
the n-th Bartnik-McKinnon soliton and the n-th non-abelian black hole solution of the 
Einstein- Yang-Mills theory for the gauge group SU(2). 
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1 Introduction 



It is by now well known that Bartnik-McKinnon (BK) solitons [|IJ and non-abelian black 
holes of the SU(2) Einstein- Yang-Mills (EYM) theory are very fragile objects. A 
linear perturbation analysis showed that these solutions are unstable and non-linear 
numerical studies M, M revealed that the instability is quite dramatic. 

The unstable modes found in j3]-|| belong to perturbations within the original BK 
ansatz. This type of instabilities will be referred to as "gravitational" , because they have 
no flat spacetime analogues. The numerical results for the few lowest members of the 
solution families led to the suspicion that the n-th EYM (soliton or black hole) solution 
has exactly n such unstable gravitational modes, a conjecture which is still not proven. 

Beside these even-parity modes there is a second class of exponentially growing 
modes, with the opposite parity, which we call "sphaleron-like" instabilities, because 
they have a similar origin as for the sphaleron solution of the Yang-Mills-Higgs system 
in flat spacetime. Heuristically, this type of instabilities is related to the existence of 
non-contractible loops in configuration space, passing through the equilibrium solution 
(at least for n odd) ||. This fact suggests that the latter are some kind of saddle points. 
For odd n's there are also paths connecting homotopically distinct vacua, differing in 
their Chern-Simons number, which pass through the equilibrium solutions [0], |§. It has 
been shown analytically that there exists at least one sphaleron-like unstable mode for 
each member of the BK family and this remains true for "generic" solitons belonging to 
general gauge groups ||, JIU| . The same has been shown for the SU (2) EYM black holes 
11 1 and "generic" non-abelian black holes for any compact gauge group |nj . Numerical 



studies [O for the few lowest members of the BK family again suggests that there are 
also exactly n sphaleron-like instabilities. It is the main purpose of the present paper to 
prove this conjecture, both for BK solitons and non-abelian black holes. 

Our demonstration is based on the study of some "dual" pairs of Schrodinger 
equations, whose partners are formally related in the same manner as in supersymmetric 
quantum mechanics. Thanks to a residual gauge group, it is possible to find the zero- 
energy solutions in sufficiently explicit form to read off the number of their nodes and 
hence, for an appropriate choice of the pair, the number of bound states, i.e., unstable 
modes. We find it quite remarkable that this can be deduced, in spite of the fact that the 
BK solutions and their black hole counterparts are not known analytically. (Rigorous 



existence proof have been given in |13|, []14|, fT5fl) 



2 Basic equations 

The basic equations for spherically symmetric SU(2) EYM fields, as well as for spherical 
perturbations of the EYM solitons and black holes, have been derived before (even for 
general gauge groups [|]). We follow the notations in Ref . [pi , where the time-dependent 
spherically-symmetric gravitational and gauge fields are parameterized as follows. 
For the metric we use Schwarzschild coordinates 

ds 2 = -NS 2 dt 2 + N~ l dr 2 + r 2 {d$ 2 + sin 2 tidtp 2 ), (1) 

and the 577(2) gauge potential is represented as 

A = aoT r dt + a\T r dr + (w — ^{r^dd — t$ sin?? dip) + w^dd + sin $ dip). (2) 
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The metric coefficients N, S and gauge amplitudes a , a x , w, w are functions of r and t, 
and T r , t&, Ty denote the spherical generators of SU (2), normalized such that [r r , t#] = r v , 
etc. 

The coupled EYM field equations, expressed in terms of the functions in (P and 
(Q), can be found in [16]. In the static case, with ao = a\ = w = 0, these equations are 
known to admit soliton and black hole solutions |jlj, characterized by the number 
of nodes of the gauge amplitude w(r). Here we repeat only the resulting perturbation 
equations for these equilibrium solutions ||, fTTf , fTB] . 



Since the parity transformation, d —> ir—fl, if — *> y?+7r, is a symmetry operation, the 
even-parity perturbations (SN, SS, Sw) and the odd-parity modes (Sao, Sax, Sw) decouple. 
The perturbation equation for the even-parity ("gravitational") modes can be brought 
into the form of a P-wave Schrodinger equation H 

- V" + (s*N 3 -^f± + 2 (§)%^ ^ (3) 

where Sw = rj(p) exp(iut), and a prime always denotes differentiation with respect to 
the radial coordinate p defined by 

^ = — (4) 

dr NS V ; 

(p(0) = for solitons and p = — oo at the horizon for black holes). 

For the odd-parity ( "sphaleron" ) modes we choose the temporal gauge, Sa = 0, 
and adopt the following notations for the perturbation amplitudes Sax, Sw: 

Sa\ = a(p) exp(iut), Sw = £(p) exp(iut). (5) 

For a and £ we obtain the coupled system, consisting of the second order equation 

/vs* 2 

- C + ~pr(w 2 - l)f + NSw'a + (NSwa)' = uo 2 ^ (6) 
and the first order equation 

^™ 2 a + ^K£-0 = ^ 2 «- (7) 

In addition, the Gauss constraint gives 

"{(^«) -2<|=0. (8) 

The following observation simplifies the further discussion of this type of perturbations 
considerably. If we multiply @ with w and take the background Yang-Mills equation 

w " = -^w(w 2 - 1) (9) 

into account, the left hand side of @ becomes a total derivative of the left hand side of 
(0), multiplied by r 2 /S. Introducing the quantities 

r 2 2 2NS 2 , x 
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we find from (||)-(|]): 

(w^x + w'Z-w?)' = u 2 w£, (11) 
w^x + w'i-wi' = u 2 X , (12) 
uj( X '-wO=0. (13) 
When uj = 0, there is a residual U(l) gauge symmetry of these equations: 

Z^Z + w n, x^x + n'h\ (14) 

where Q is a function of p. This will be of crucial importance. 

3 The number of odd-parity sphaleron instabilities 

In this section we prove that the number of negative (u 2 < 0) modes of the system 



(11)- (13) for the n-th BK solution, as well as for the n-th EYM black hole, is exactly n. 
We call this type of instabilities sphaleron-like, because of their similarity with the flat 
spacetime Yang-Mills-Higgs sphaleron instabilities [0. 

For uj 0, Eq . (|TT|) is a consequence of (|12"D and (|T3l) , i.e., 

<-tD'(=(wV-w 2 )x, (15) 
X ' = w{. (16) 



As a technical tool we keep equation (jig ) also for uj = and call it strong Gauss 
constraint; it plays the role of a gauge fixing condition. It suffices to determine the 
number of uj 2 < modes of the system ([15]), (|T6|). (The role of the Gauss constraint is 
discussed more extensively in fllPl ). 

We begin our analysis by eliminating £, with the result 

- 0" + {]p\w 2 + 1) + 2 ^ 2 j = u 2 ^ (17) 

where = x/ w - F° r a smooth function w without zeros this would be a regular 
Schrodinger equation with a purely positive spectrum {u 2 > 0). This suggests that 
the number of unstable modes of the n-th equilibrium solution is indeed related to the 
number n of nodes of w. 

Since the differential equation fllTD has singularities at the node positions of w, we 
try to pass to a "dual" regular Schrodinger equation, for which the number of bound 
states can be read off. How to achieve this is suggested by the following considerations. 

First, we note that Eq.([T7|) has obvious u = solutions which are just pure gauge 
modes (see (|14D ), 

0o = n'/w 7 2 , (is) 



satisfying the strong Gauss constraint (|16| ) with £ = wfl: 

/MX' 



, , = w 2 n. (19) 

T 
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Any such solution allows us to factorize the differential operator in (|17|): 



,/\ 2 \ ^2 



^ + i 7>2 + 1) + 2 ^ = _£_ + £ = 0+0 -, (20) 



with 



Q* = t4- — rr- (2D 



dp 



Using (|19D we find 



where 



Q ± = Ty" + — + w 2 Z, (22) 

Z = -QA, A = 7 2 /fi'. (23) 
For later use we note that (fL9|) is equivalent to 

A/A = w 2 Z, (24) 

and thus Z satisfies the nonlinear differential equation 

Z' = w 2 Z 2 - 7 2 . (25) 

With a standard reduction, the most general solution of the second order linear 
equation (|T9|) is 

Q = c 2 A + tt (cx + c 2 £ w 2 A 2 dp^j , (26) 

where Q is a special solution, A = 7 2 /0', and ci, c 2 are real constants. This gives 
immediately for Z in Q2"B] ) 

A 2 

Z = -OA p — : . (27) 

Ci/c 2 + / w 2 A 2 dp 
Jo 

Since the differential equation ([D]) is, for any Q in fl26|) , identical to 

Q+Q-0 = w 2 0, (28) 

we are invited to pass from to 

i) = Q~<J). (29) 

For w ^ this has the unique inverse 

= \Q + ^ (30) 
and by applying Q~ we obtain the "dual" eigenvalue equation 

Q-Q + ij = u 2 i>. (31) 
The differential operator on the left is 

Q-Q+ = ~^ + U, (32) 
4 



with 

U = -7 2 (3w 2 -l) + 2(w 2 z)'. (33) 

In this potential the nodes of w are no more dangerous and we have the chance that fl3~l|) 
becomes a "regular" Schrodinger equation if Q in ( PU] ) is chosen appropriately. 

Before addressing this problem, we determine the zero-energy solutions of ([31~D. 
These satisfy Q + ipo = 0, i.e., 



4>o w 



- w + W *Z = ^, (34) 



thus 

1 A 

ipo = — = wA = w jp . (35) 

ci + c 2 / w 2 A 2 dp 
Jo 

An alternative form is obtained from (24), (35): 

ipo = wexp ^ w 2 Zdp s j . (36) 

This looks good, because the factor w would count the number of nodes of this zero- 
energy solution if the second factor on the right hand side could be chosen to be every- 
where regular. 

A detailed discussion of the function Z in (^7|), whereby fl is a special solution of 
the linear equation (fL9|), will be given elsewhere. Here we only state that one can choose 
Q such that 

~ coin nt ~ \ \ 

Q= ^1 + 0(1) asr — 0, Q= 4 + 0(4) as r ^ oo (37) 

for the BK solitons, and 

(l = x + 0(x 2 ) as x — > 0, (l = const x + 0(1) as x — > oo (38) 

for the EYM black holes, where x = (r — r^/r^, being the Schwarzschild coordinate 
of the horizon. 

For the regular solutions, one obtains from fl3?D, (|27]) (with an appropriate constant 



crijc\ in (27|)) that Z exists everywhere on the half-line with the following behavior near 
p = 0, oo: 

Z = - + 0(p) asp->0, Z = -4 + 0(^^) as p -> oo. (39) 
p p 2 p^ 

This assures that the potential [/ in ( |33|) is everywhere bounded. Indeed, the first term 
for [/ behaves as 2/p 2 for p — * 0, and this singularity is exactly canceled by the second 
term in U involving the function Z. 

In summary, the "dual" equation ([H]) can be chosen such that it becomes an .S-wave 
Schrodinger equation with an everywhere bounded potential. In addition ( |39| ) assures 
the following behavior for the zero-energy solution ipo: 

ipQ ~ p as p — > 0, ip$ ~ — as p — > oo, (40) 

P 
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which implies that ipo G L 2 [0, oo). A well-known theorem ensures then that the number 
of nodes of ipo, which is n for the n-th BK soliton solution, is also equal the the number 
of bound states of ([H]), and thus equal to the number of sphaleron-like instabilities. 



Similarly, for black holes, the special solution (|38l) , together with an appropriate 
choice of c^jcx in (0), ensure that Z has the following behavior near p = ±00: 

' 0{\) as Z -> -^ + 0(^^) as p -> 00. (41) 



f 2 (r h )p p 2 ' p 2 p 3 

This yields 

1 1 

ijj ~ - as p — > —00, ^0 ~ as P - >■ °°? (42) 
P P 2 

thus ^0 G L 2 (— 00, 00), implying that for the n-th EYM black hole (independent of its 

event horizon size) there exist precisely n unstable odd-parity modes. 



4 Concluding remarks 

Looking back at the proof, the following elements have been crucial. Thanks to the 
residual gauge symmetry (III) we know zero energy solutions of the original differential 
equation fll7f) and we can thus factorize the differential operator. We can choose this 
factorization such that the "supersymmetric partner" is a well-behaved Schrodinger equa- 
tion which has a normalizable zero-energy solution. For this solution, we can determine 
the number of zeros and hence read off the number of bound states. 

Unfortunately, because 77 in Eq.(|3|) remains invariant under the residual gauge 
transformation, this method does not work for the even-parity gravitational instabili- 
ties, whose number is presumably the same. 
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